Random Isotropic Structures and Possible Glass Transitions 
in Diblock Copolymer Melts 
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We study the microstructural glass transhion in diblock-copolymer mehs using a thermodynamic 
repUca approach. Our approach performs an expansion in terms of the natural smallness parameter- 
the inverse of the scaled degree of polymerization iV-which allows us to systematically study the 
approach to mean-field behavior as the degree of polymerization increases. We find that in the 
limit of infinite chain length, both the onset of glassiness and the vitrification transition (Kauzmann 
temperature) collapse to the mean-field spinodal, suggesting that the spinodal can be regarded as 
the mean-field signature for glass transitions in this class of microphase-separating system. We also 
study the order- disorder transition (ODT) within the same theoretical framework; in particular, we 
include the leading-order fiuctuation corrections due to the cubic interaction in the coarse-grained 
Hamiltonian, which has been ignored in previous studies of the ODT in block copolymers. We find 
that the cubic term stabilizes both the ordered (body-centered-cubic) phase and the glassy state 
relative to the disordered phase. In melts of symmetric copolymers the glass transition always occurs 
after the order-disorder transition (below the ODT temperature), but for asymmetric copolymers, 
it is possible for the glass transition to precede the ordering transition. 

PACS numbers: 61.41. -|-e, 64.70.Pf, 64.60.My, 64.60.Cn 
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I. INTRODUCTION 

Block copolymers are macromolecules built with blocks 
of chemically distinct monomers. Melts of block copoly- 
mers are attractive from both theoretical and experimen- 
tal standpoints, as they undergo microphase transitions 
and produce diverse ordered microstructures 0, S S 13 • 

The simplest block copolymer is the AB diblock 
copolymer made of two types of monomers A and B. 
Below the order-disorder transition (ODT) temperature, 
a diblock-copolymer melt can exhibit rich mesophases 
[3, including body-centered-cubic (bcc), hexagonally or- 
dered cylinder (hex), lamellar (lam), and several bicon- 
tinuous (e.g., gyroid) structures. Experimentally these 
structures have been identified using transmission elec- 
tron microscopy (TE]\^ Q , small-angle neutron scat- 
tering (SANS) [aja, and dynamic mechanical mea- 
surements 1^ I9I llOf. Theoretically these structures are 
well described by the self-consistent-mean-field theory 

iniiii. 

Generally these periodically ordered structures are ex- 
pected to be the thermodynamically equilibrium states 
. However, they are difficult to attain either in exper- 
iments or in computer simulations Bates and 
co-workers 0, Q found that quenching a nearly symmet- 
ric diblock-copolymer melt without symmetry-breaking 
external field, such as reciprocal shearing, generally re- 
sults in isotropic, locally microphase-separated structures 
with a characteristic length scale of the radius of gyra- 
tion of the polymer. In addition, such structures were 
also obtained (as a rule) in dynamic-density-functional 
calculations Ulllji,]. The ordering kinetics in these 
random structures are very slow, suggesting that they are 
metastable states corresponding to free-energy minima. 



It is therefore quite possible that the ordered phases, 
though energetically favored, are not easily reached due 
to the kinetic trapping caused by the presence of a 
large number of metastable free-energy minima. These 
metastable states correspond to the locally microphase- 
separated states without long-range order. 

Dynamic mechanical measurements by Bates and co- 
workers on melts of both symmetric, lamellae-forming, 
and asymmetric, hex-forming PE-PEP [partially deuter- 
ated poly(ethylenej)ropylene)-poly(ethylethylene)] 
copolymers revealed 0, 13 that the system may be frozen 
in random structures upon a deep quench. Comparing 
the quenched sample with the slowly-supercooled sample 
and the shear-ordered sample, they found that the 
quenched sample exhibits very slow relaxations and 
extraordinarily large elastic moduli at low frequencies; 
but the supercooled sample behaves more like the dis- 
ordered melt continuously extended to below the ODT 
temperature. Balsara and co-workers 0, studied 
the grain structure of asymmetric, hex-forming PI-PS 
(polyisoprene-polystyrene) melt by light scattering, 
SANS, and rheological measurements. Similar to the 
findings of Bates et ai, they found that upon a deep 
quench, randomly microphase-separated structures are 
obtained, which do not appear to evolve towards the 
equilibrium structure with long range order within the 
time scales of the experiments. Besides these, Pochan 
et al. jl7j found randomly oriented wormlike cylinder 
structures in an I2S [polyisoprene(I)-polystyrene(S)] star 
copolymer system,. 

The above results suggest that the ordering process in 
block-copolymer melts follow a two-step mechanism: a 
fast step in which unlike monomers locally phase sep- 
arate into random, macroscopically isotropic structures 
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with domains of the size of a single polymer, followed 
by a domain coarsening (or growth) step in which local 
defects in the random microstructures anhilate and long 
range order is developed. The second step is generally 
much slower than the first and most likely involves ac- 
tivated processes. Therefore a rapid deep quench can 
result in randomly microphase-separated structures that 
are kinetically trapped and unable to develop long range 
order within normal laboratory time scales. 

This two-step mechanism is in fact consistent with 
the thermodynamic two-step scenario implicit in the 
Fredrickson-Helfand (FH) fluctuation theory for diblock- 
copolymer melts |[18j (which only applies to symmet- 
ric or nearly symmetric copolymers). Instead of the 
featureless background as assumed in the random-ficld- 
approximated structure factor of Leibler |0| , the FH the- 
ory suggests that when the temperature approaches the 
ODT, the disordered state is a fluctuating, heterogeneous 
structure consisting of locally A- and _B-rich domains, 
which then orders into periodic mesophases upon further 
cooling. 

The most dramatic manifestation of the first step 
is the existence of disordered-spherical-micelle state in 
highly asymmetric copolymer melts, which almost has 
the appearance of a distinct phase between the fea- 
tureless disordered phase and the bcc-ordered phase 
[2Sp,|2l|2l|2l|2l|2iiil. The micelle state was first 
predicted by Semenov [2g. Recently Dormidontova and 
Lodge |29| | extended the Semenov theory by including 
the translational entropy of the disordered spherical mi- 
celles and predicted a phase diagram that is in qualitative 
agreement with experiments. More recently, Wang and 
co-workers p30| examined the nature of the disordered 
spherical micelles and their connection to concentration 
fluctuations using the self-consistent-field theory. Tak- 
ing a nucleation perspective, these authors showed that 
the disordered micelles are large, localized concentration 
fluctuations through a thermally activated process. 

In this work, we study the metastable states consist- 
ing of random structures in block-copolymer melts and 
address the possibility of glass transition using a ther- 
modynamic replica approach. This approach was first 
proposed by Monasson "31] and subsequently employed 
by a number of authors in studying structural glass tran- 
sitions [33, 113,0, HE 113. In this framework, the onset of 
glassiness is identified with broken ergodicity , which 
occurs as a result of the appearance of an exponentially 
large number of metastable free-energy minima |3Sj . The 
broken ergodicity is manifested through a nonvanishing 
long-time correlation (here manifested as the cross replica 
correlation function), whose first appearance defines the 
onset temperature of glassiness Ta (also called the dy- 
namic glass transition temperature |3ll l35l|). An equiv- 
alent Kauzmann temperature Tr- as in molecular liquids 
|39] can also be defined as signaling the complete vitrifi- 
cation of the random structures. 

The possibility of glass transitions in bicontinu- 
ous microemulsions-a system closely related to diblock 



copolymers-was recently examined by Wu, Westfahl, 
Schmalian and Wolynes using both a dynamic 

mode-coupling theory and the thermodynamic replica 
approach. There authors have also studied glass transi- 
tions in the Coulomb-frustrated-magnet model using the 
replica method with a self-consistent-screening approx- 
imation '3^, '3^ and, more recently, a local-field calcu- 
lation 36]. Both the microemulsion and the Coulomb- 
frustrated-magnet systems belong to the general class of 
models first proposed by Brazovskii [4]| , featuring the ex- 
istence of low-energy excitations around some finite wave 
number and the formation of microphase-separated 
structures with length scales ~ l/q-m at low tempera- 
tures. These studies showed that as a result of the large 
degeneracy in ground states (s^, a glass transition can 
occur when the ratio of the correlation length of the 
system to the modulation length 2T:/qm exceeds some 
critical value. Similar conclusions were also obtained by 
Grousson et al. using the mode-coupling theory. 



Our work follows a similar approach to that employed 
by Schmalian and co-workers 34, 35]. However, we per- 
form calculations specifically for the block-copolymer sys- 
tem by taking advantage of the natural smallness param- 
eter (the inverse of the scaled degree of polymerization, 
N); this allows us to study how the glass transitions are 
affected by increasing the chain length of the polymer 
when the system gradually approaches the mean-field 
limit. An important conclusion of our work is that in the 
limit of infinitely long chains, both the onset of glassi- 
ness and the Kauzmann temperature coincide with the 
mean-field spinodal of the disordered phase. Therefore 
the spinodal is the mean-field signature for the glass tran- 
sition in the block copolymer system; the same conclu- 
sion is likely to hold in general for microphase-separating 
systems. Another feature of our work is the inclusion 
of the order-disorder transition in the phase diagram. 
This is important because it places the glass transition 
in proper relationship to the ordering transition. We find 
that, for symmetric, lam-forming copolymers, the glass- 
transition temperatures are below the ODT temperature, 
while for asymmetric, sphere-forming copolymers, the 
onset of glassiness can precede the ODT into the bcc 
phase. On a technical point, we propose a method for 
incorporating fluctuations due to the cubic interaction 
in the Brazovskii model, using a renormalization scheme 
motivated by the 1/n expansion of the n- vector model 
in critical phenomena. The effects of these fluctuations 
have not been addressed in any of the previous studies 
0,000,0] on block-copolymer systems. We find 
that in the leading-order approximation these fiuctua- 
tions stabilize both the bcc phase and the glassy state. 
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II. MODEL AND SOLUTION 



where 



A. Model description 

We consider the melt of AB diblock copolymers of de- 
gree of polymerization N = Na + Nb and block com- 
position / — Na/N. The monomer volume v and Kuhn 
length b are taken to be equal for both monomers. We 
describe the thermodynamics of the system using the 
random-field-approximated (RPA) free-energy functional 
with local approximations for the cubic and quartic in- 
teractions |l8i llfll li^ as the Hamiltonian 

^ ^' ~q)l2{q,~q)(l){q) 
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where the order parameter = pa{x)v — / is the density 
deviation from the mean value. Throughout the paper 
we take fc^T = 1 except for our discussion of the ther- 
modynamic approach to the glass transition in Sec. Ill (Jl 
To simplify the notation, we use plain letters (x, q, etc.) 
to denote position and wave vectors; when the plain let- 
ter is used to denote the magnitude of the wave vector 
(wave number), the context should make it clear. 

Near the mean-field spinodal 72(9, —q) can be approx- 
imated as 



l2iq,-q) = j{q^Nb^ 

where X-ZV is the Flory-Huggins interaction parameter be- 
tween A and B blocks, (x^)s is its value at the spinodal, 
and c is a parameter independent of N. (x^)s^ c, and 
qm are functions of / and N, which can be calculated 
using the RPA theory of Leibler [l^. Note that Eq. (HJ 
as a Hamiltonian is applicable to a broad class of copoly- 
mer systems, including multiblock copolymers |43j| and 
copolymer/homopolymer blends p[3 |. where the depen- 
dence on chain architectures, block compositions, and 
volume fractions of copolymers can be incorporated into 
the parameters (x-^)sj qm, etc. Therefore our results on 
diblock-copolymer systems should be qualitatively appli- 
cable to these systems as well. 

The degree of polymerization, N, plays the role of 
Ginzburg parameter, which controls the magnitude of 
fluctuations To highlight this feature, we nondi- 

mensionalize the lengths and wave numbers by the ideal 
end-to-end distance of the polymer: x = x/{^/Nb), q = 
qVNb, qm = qmVNb, and concurrently rescale the order 
parameter as (j}{x) = (j){x)cqm, (t>{q) = 4'{q)cqm/ {yfNbY . 
Now the Hamiltonian [Eq. (^] becomes 
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(4) 
(5) 
(6) 



(2) 



The scaled couplings rj and A are, respectively, the same 
as NT3 and A^r4 defined in Ref. 0|. For notational 
simplicity, we drop the overbars on the variables and the 
order parameter henceforth. 

We point out that, although the parameters in Eq. (0 
are written in molecular terms, this model is best in- 
terpreted as phenomenological. The random-phase ap- 
proximation used in deriving the Hamiltonian, the ap- 
proximation of higher-order interactions as spatially lo- 
cal, and the truncation at quartic order in the order- 
parameter expansion-all introduce inaccuracies whose ef- 
fects are difficult to evaluate In particular, the 
order-parameter expansion to quartic order is not jus- 
tified for strongly asymmetric block compositions as 
chain stretching effects become important and the weak- 
segregation assumption no longer holds "s^. However, 
we note that taking Eq. Q as the Hamiltonian, one can 
reproduce the experimental phase diagram of microphase 
transitions qualitatively at all compositions, including 
the disordered spherical-micelle states at very asymmet- 
ric compositions, as the state-of-the-art self-consistent- 
field theory. Therefore while the quantitative accuracy 
of our theory may not be reliable, we expect that most 
of our predictions should be qualitatively correct. Such 
an expectation is further boosted by the general success 
of the Fredrickson-Helfand theory [also using Eq. |5J as 
the Hamiltonian] in capturing many key features of the 
physics of diblock-copolymer melts at length scales com- 
parable to or larger than the size of the polymer chain. 

In addition, studying glass transitions in the system 
described by Eq. ||2J) is of intrinsic theoretical value, as 
Eq. ^ corresponds to the weak-coupfing fimit of the 
Brazovskii model. Therefore our results elucidate the 
physics of systems in the Brazovskii class in this limit. 

Finally we notice that the parameter iV^/^ = iV^/^6^/u 
(henceforth referred to as the "chain length" ) is a natural 
combination emerging in any study of the fluctuation ef- 
fects in polymer melts, which gives the number of other 
chains within the spatial extension of a single polymer 
chain [TsL l5ll |. N plays a role similar to 1/h in quan- 
tum field theory [S^l-controlling the magnitude of fluc- 
tuations. In the limit of iV ^ 00, mean- field behavior is 
recovered. For systems with large but finite N we can 
apply a systematic loop expansion using l/N^^^ as the 
smallness parameter. 

The presence of the 7V^/^ factor in front of the Hamilto- 
nian also has important consequences on the free-energy 
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barriers separating the multiplicity of free-energy min- 
ima. In the mean-field approximation, we expect that 
the free-energy barriers should be proportional to this 
factor. For long polymer chains, the barriers can be 
much larger than the thermal energy, resulting in slow 
relaxations from these free energy minima to the lower- 
free energy ordered phases and between the metastable 
minima themselves. This justifies the application of the 
energy-landscape theory of glass transitions in polymer 
systems. 

B. Ordered states and order-disorder transition 

Our current understanding of the effects of fluctuations 
on the ODT in block-copolymer melts is largely based on 
the Brazovskii-Leibler-Fredrickson-Helfand (BLFH) the- 
ory 0, l4ll | . This theory uses the self-consistent Bra- 
zovskii approximation (a Hartree-type approximation) 
for the quartic interaction and ignores fluctuations due 
to the cubic interaction. Therefore, strictly speaking, it 
is only valid for symmetric or nearly symmetric block 
copolymers where cubic interaction is small (see our dis- 
cussions at the end of this subsection). Here we extend 



this theory to include the leading-order one-loop correc- 
tion from the cubic interaction, which accounts for the 
fluctuation effects due to asymmetry in the copolymer 
composition. This improved theory should give more ac- 
curate predictions on the ODT in asymmetric copolymer 
melts (and other asymmetric systems) and, more impor- 
tant, enables a consistent comparison with the glass tran- 
sition in the same system, where the cubic term is shown 
to play a dominant role. 

As in previous weak-segregation theories 0, 0] , we 
adopt the single-mode approximation for the periodic mi- 
crophases, representing the density wave by 

(p{x) = [exp{iQj ■ x) + exp{—iQj ■ x)] , (7) 
j 

where a is the magnitude of the density wave and Qj{l < 
j < n) are the first set of vectors on the reciprocal lat- 
tice of the periodic structure of the ordered microphases 
[TsL Hoi . Now we introduce the fluctuation field around 
the minimum, il'{x) = 4'{x) — (p{x), and perform an ex- 
pansion of the Hamiltonian H[(l)] in Eq. (jSJ around ip. 
The fluctuation part of H[(j>] is 



2 J (27r)3"' ' 3! J "---^"•^ ■ 4! 



(27r)9 

^ j d pid P2d^p3d PA [4^(-p^)^(-p2)^/,(p3)^(p4) + QLp{pi)ip{p2)lp{p3)lp{pi)] 5^(pi +P2 +PZ+Pi). (8) 

The linear term of "0 vanishes because ip is at the minimum of the Hamiltonian. For the quadratic term we only keep 
the dominant isotropic part 

D{qr^^g{qr^+n\a\ (9) 

which is defined as the shifted bare propagator. 

The free energy (effective potential) of the microphase-separated system is given by 



F[ip\ = ~N-^/'^\n(eyiY>^-N^/'^H[(j)\Y^ 

^ H[ip]-N-^''^\n{e:^Y'[-N^'^AH[il}-,ip\yj. (10) 

Here the free energy is scaled by 7V~^/^ such that the mean-field part H[lp\ is independent of TV and reduces to 
the Leibler free energy The second term in Eq. ((Tn|l contains corrections due to the fluctuation part of the 

Hamiltonian [Eq. JHJ] . In the one- loop approximation we have 



^M=^M + ^Trlng--A 



■ Gniq) 



(27r)3 



UN J (27r)6 



where Gniq) is the Hartree-renormalized propagator determined from 



A f d^k 

5Hiq)-^-Diq)-^ + ^l-—g,ik). (12) 



5 



Our one-loop approximation is slightly different from the conventional diagrammatic expansion; details are discussed 
in Appendix IBl 

Under this approximation the renormalized correlation function is given by 



\dip{q)6ip{-q) 

2Ari/2 J (27r)3^ 



n2 I- j3l 

Sniq)-' - / 77;Z^gH{k)gH{q k). (13) 



In the replica calculation Q{q) gives the renormalized di- 
agonal correlation function in the replica space. 

The second term in Eq. H13|l . corresponding to the one- 
loop cubic diagram, was absent in previous theories on 
the ODT, as it was shown to be subdominant to the 
Hartree term [the first term in Eq. near the mean- 
field ODT for asymmetric copolvmers |4lLl5^ . The argu- 
ments for ignoring this term no longer hold for the super- 
cooled disordered phase (below the ODT temperature), 
and here we need a free-energy function that remains 
valid even below the mean-field spinodal temperature. 
Therefore the one-loop cubic term cannot be dropped as 
in Refs. and |0. Also in the 1/N expansion em- 
ployed here (equivalent to a loop expansion), the one- 
loop cubic term is of the same order as the Hartree term 
and their numerical values are comparable in the part 
of the phase diagram of interest, except for nearly sym- 
metric compositions when the cubic term is small |74j . 
Furthermore, as shown in Appendix fXlthe corresponding 
term is the leading term in the self-consistent equation for 
the cross-replica correlation function. Earlier work also 
showed that it is the leading term that generates long- 
time correlations in the mode-coupling theory for glass 
transitions js^ . We therefore include the one-loop cubic 
term in our treatment of the ODT to have a consistent 
comparison with glass transition. 

C. Random structures and glass transition 

Traditionally two different approaches have been devel- 
oped to study frustrated systems with quenched disorder. 
The dynamic approach, most notably the mode-coupling 
theory |5^], focuses on dynamic correlation functions 
(e.g., the Edwards- Anderson order parameter defined as 
the long-time spin-spin correlation function in the Ising- 
spin-glass model l56f) and characterizes the glassy state 
with nonvanishing long-time correlations and broken er- 
godicity. On the other hand, the equilibrium thermo- 
dynamic approach, including the density functional ap- 
proach 57] and the replica approach describes glass 
transitions in terms of the energy-landscape features of 
the system [s^. The connection between these two ap- 
proaches was explicitly demonstrated in the mean-field 



spin-glass models |58|, |59j where it was shown that these 
two approaches yield consistent predictions. We now 
briefly describe the essential concepts in the thermody- 
namic approach. 

The central assumption in the thermodynamic ap- 
proach is that the dynamic behavior of glass- forming sys- 
tems reflects the un derly ing free-energy-landscape fea- 
tures [sil Is^ Is^ l60l l6l| . At high temperatures, there 
is only one minimum corresponding to the uniform liq- 
uid state. As temperature decreases, multiple metastable 
minima begin to appear that are separated by sizable 
activation barriers, and below some temperature T^, 
the number of these minima becomes thermodynamically 
large, giving a finite contribution to the partition sum of 
these "disjoint" metastable states and generating exten- 
sive configurational complexity manifested in a nonvan- 
ishing configurational entropy ||61| . This signals the onset 
of glassiness or broken ergodicity js^j in the sense that 
within times scales of typical liquid relaxations, the sys- 
tem is trapped in these metastable free-energy minima; 
transitions between the minima, however, can still occur 
through activated processes |59|. Dynamically, one ex- 
pects a significant slowing down of structural relaxations, 
often accompanied by the appearance of long plateaus in 
the time correlation functions [39l Is^ . Complete vitri- 
fication occurs at a lower temperature Tr-, below which 
the system is dominated by one or less than an expo- 
nentially large number of deep free energy minima; ther- 
modynamically, this is signaled by the vanishing of the 
configurational entropy. Tk is often termed the ideal 
glass transition temperature and is conceptually identi- 
fied with the underlying thermodynamic glass transition 
at which the viscosity of the supercooled liquid diverges 
Hi, 39, 60, 62]. 

Recently Monasson (3l| proposed a replica method 
which allows explicit implementation of the thermody- 
namic approach for studying structural glasses result- 
ing from self-generated randomness. Using this method 
Westfahl et al. [l^ successfully predicted the glass tran- 
sitions in the Coulomb-frustrated-magnet model. Here 
we adopt this approach to study the glass transition in 
block-copolymer melts. 

Following |3lj] , we introduce an external pinning field 
and calculate the pinned free energy of the system, F[C] 
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F[C] = -kBTlnZlC] = -fcsTln / X>0exp 



(14) 



where a > is the couphng between the pinmng field and the order parameter. [In Eqs. p4|l -H2(J () . we reintroduce the 
ksT factor in order to allow explicit temperature derivatives.] The effect of C is to locate the basins on the free-energy 
landscape. The coupling constant will be taken to be infinitesimally small at the end and serves as a convenient device 
for breaking ergodicity-localizing the system into separate basins. Its role is similar to that of the infinitesimal field 
that breaks the up-down symmetry of the Ising model below the critical temperature. One can show that the minima 
of F[Q coincide with those of the effective potential of 7i[0] as a ^ 0; proof is given in Appendix IdI Thus ( serves 
as a running index for labeling different basins on the free-energy landscape, and sampling the configuration space 
of ^ gives information on the metastable free-energy minima (the energy minima with their location fiuctuations) of 
the system. Therefore one can use F{() as an "effective Hamiltonian" for the metastable free energy minima and 
compute the "quenched-average" free energy 



F 



JVCF[C]e^p{-F[C]/kBT} 
JVCcxp{-F[C]/kBT} ■ 



(15) 



If the system is fully ergodic, one can verify that F is equal to the equilibrium free energy 



F = -kBT\n j V(j)eyi^{~'H[(l)\/kBT} 



in the thermodynamic limit as a 
Their difference 



0+. However, when ergodicity is broken limQ^o+ P can be different from F. 



F -F ^TSc 



(16) 



defines the configurational entro py that measures the configurational complexity due to an exponentionally large 
number of metastable states |3ll l35l IstI l58l |. In the thermodynamic approach, Sc jumps discontinuously from zero 
to an extensive finite value at Ta, implying broken ergodicity due to disjoint metastable states; Sc decreases upon 
further cooling and vanishes at T^, when the system becomes completely vitrified. 
To calculate Sc it is convenient to introduce the "replicated" free energy 



Fm. = - lim 



kBT 



knT 



■ In Zr, 



(17) 



where T/m is introduced as the effective temperature conjugate to F[Q. F and Sc are obtained from Eq. 
straighforwardly as 



F = 

Sc = 



d{mF„ 



dm 

dFrr 



d(T/m) 



1 dF„ 



m—1 



T dm 



(18) 
(19) 



m— 1 



When m is an integer, Zm in Eq. 1)17(1 can be simplified by introducing m copies of (j) and integrating out the C field, 
which gives 



~ lim 



Vcpa CXp < 



a— 1 



2mkBT 



m „ 



(20) 



l<a<fc<m 



where a, 6 are replica indices. Equation H20() has the same form as the replicated partition function for a random 
system with quenched disorder 58], although here we are interested in the physical limit corresponding to to = 1. 

To characterize the physical states of the system, we introduce the (renormalizcd) correlation functions Q(ci) = 
{4>a{<l)4>a{—<l)) and ^(g) = {4>a{<l)4>b{—<l)) a=/^b- ^il) is the normal physical correlation function of the system, whereas 
^{q) measures the correlation between different replicas. It has been shown that f^{q) is equivalent to the long-time 
correlation function in the conventional mode-coupling approach [sil l35l l54j . At high temperatures the system is 
ergodic, and in the limit a ^ 0'*' different replicas are not coupled; thus, J-{q) = 0. When ergodicity is spontaneously 
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broken, different replicas become coupled even in the limit a ^ 0+, and J^{q) 7^ 0. Using T{q) as the order parameter 
for ergodicity breaking, we can define the onset of glassiness Ta as the temperature when there first appears a solution 
with !F{q) 7^ 0. Ta defined in this way coincides with the dynamic-transition tem per ature in mean-field-spin-glass 
models characterized by the appearance of drastically slow dynamic relaxation [sil iBSj . 

To obtain the replica free energy defined by Eq. (|17|l . we adopt the self-energy approach [s^ and express the effective 
potential Fm as a functional of bare and renormalized correlation functions 

F„, [G] = - (ilrlnG-i + ^Tr (D-^G) - r2[G]l . (21) 
m I 2 2 I 



Here D and G are bare and renormalized correlation 
hmctions, with 

G = {g-T)l + T-E, (22) 

where Eab — 1 and lab ~ Sab- [Henceforth we use bold 
face uppercase letters (G, D, S) for the matrices of func- 
tions in the replica space, plain uppercase letters with 
subscript indices (Gab, Dab, etc.) for the matrix ele- 
ments, g and J- are reserved for the renormalized di- 
agonal (physical) and cross-replica correlation functions 
respectively] 

Dab is the replica-symmetric bare correlation function, 
Dab{q) = g{q)Sab, with g{q) given from Eq. ® 

The self-energy functions Y^ab arc defined by 

S = G"^-D"i (23) 
and obtained through variation of Fm'- 



r2[G] contains all two-particle-irreducible (2PI) dia- 
grams, which is evaluated perturbatively. Detailed cal- 
culations are given in Appendix IXI 

Taking the inverse of G defined in Eq. I|22|) . we find 
that the self-energy from Eq. H23fl takes the form 



Safc = (S^ - T,jr)5ab + SjF, (25) 



where 



S5(g) = -.9(9)-', (26) 



Assuming that the momentum dependence of self- 
energy functions '^g{q) and Tij^{q) is negligible compared 
with g{q), we can approximate the renormalized diagonal 
correlation function as 



a™ e(Q)-^ « i ( J - 1) +ro + Se(g,„)9,-2^i(J-l) + r. (28) 
And the off-diagonal correlation function J-{q) takes the form 




Equations. H23() and H24|l give the self-consistent equations for Q and T (algebraic equations for r and s in our 
case). Solving these equations we obtain a normal replica-symmetric solution with r = s and a replica-symmetry- 
broken solution with r < s below the dynamic-transition temperature Ta [corresponding to some {xN)^ in our 
diblock-copolymer model] . 
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The configurational entropy is obtained from Eqs. (|19() and 121() to be 



(2^ 



In 



Giq)) 



_d_ 

dm 



(32) 



One indeed finds that Sc becomes extensive below Ta 
and decreases to zero at T = Tk < Ta', Tk determines 
the Kauzmann temperature or the thermodynamic glass 
transition defined above. 



III. RESULTS AND DISCUSSION 
A. Glass transition 

Glass transitions in the Coulomb-frustrated-magnet 
model have been addressed by several groups in recent 
years 0^1^1138, 42, 63, 64, 65, 66]. These studies 
establish that in this model glass transitions are possi- 
ble and could be kinetically favored. However, all these 
studies focus on the strong-coupling regime, and except 
in Ref. '36], the asymmetric cubic interaction has been 
ignored. The block-copolymer system we are studying 
belongs to the same universality class as the Coulomb- 
frustrated-magnet-both are examples of the Brazovskii 
model. But for long chains our system corresponds to 
the weak-coupling regime of the Brazovskii model. Fur- 
thermore, the presence of the cubic interaction, reflecting 
compositional asymmetry in the copolymer, is the rule 
rather than exception. It has a strong effect on the ODT 
and the glass transition, as we will discuss in this work. 

We start with the glass transition. Figure ^ shows 
the transition lines for two chain lengths N = 10^ and 
iV = 5 X 10^. The dotted line represents the mean- field 
spinodal; dashed lines represent the Kauzmann tempera- 
ture (or the thermodynamic glass transition temperature 
|35|') Tk- The dynamic glass transition temperature Ta 
is found to be close to the Kauzmann temperature Tk 
on the scale of this figure in both cases, so we do not 
present Ta here and only include it in Figs. 2 and 3. In 
the energy-landscape theory of glass transitions, Ta sig- 
nals the onset of glassy behavior (e.g., slow dynamics), 
whereas Tk represents the limit of supercooling below 
which the system becomes vitrified 3^. (Note that we 
use the term "temperature" even though the phase dia- 
gram is presented in terms of the Flory-Huggins interac- 
tion parameter x^'^ the actual temperature can be deter- 
mined from the temperature dependence of ■) These 
results show that in diblock-copolymer melts, glass tran- 
sitions occur at finite temperatures at any chain compo- 
sition /. The narrow gap between Ta and Tk suggests 
that the system becomes vitrified right after the onset of 
glassiness. Furthermore, as the chain length increases, 
both Ta and Tk transitions approach the mean-field 
spinodal. This latter result is consistent with our antici- 
pation that a large number of inhomogeneous metastable 



free-energy minima emerge as the system approaches the 
mean-field spinodal. 

Figure^also shows the full crossover from nearly sym- 
metric copolymer, whose glass transitions are dominated 
by the quartic coupling, to highly asymmetric copolymer 
dominated by the cubic coupling. (We again remind the 
reader that the results for highly asymmetric block com- 
positions should only be taken as qualitatively but not 
quantitatively valid.) For symmetric or nearly symmetric 
copolymer, it is well known [TsL l6^ that the mean-field 
spinodal is destroyed by fluctuations and the disordered 
phase is always locally stable. Also the transition from 
the disordered phase to the lam phase is a first-order 
transition with rather complicated (and probably slow) 
kinetics (46- 68 J . Therefore a deep quench without an- 
nealing can result in the trapping of the system in ran- 
domly microphase-separated structures; these structures 
represent the glassy state captured here. This scenario 
is consistent with the experimental observations of Bates 
et al. 0, where they studied the mechanical properties 
of three different samples: a rapidly quenched sample, 
a slowly supercooled sample, and a shear-oriented sam- 
ple. By analogy to molecular liquids, these three samples 
can be likened to the glassy state, the supercooled-liquid 
state, and the ordered crystalline state, respectively. The 
quenched sample in this study exhibits solidlike responses 
at low frequencies while the supercooled sample has typ- 



XN 



Kauzmann 
Spinodal 



N = 10'' 



= 5 X 10"^// 
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FIG. 1: Glass transitions in diblock-copolymer melt. Dashed 
lines are for the Kauzmann temperature and the dotted line 
for the mean-field spinodal. The upper dashed line is for chain 
length iV = lO"*; the lower one for iV = 5 x 10''. Since Ta and 
Tk are very close, only the Tk transition is shown here. 
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ical liquidlike responses. 

We notice that in going from symmetric to asymmetric 
compositions on either side, the transition lines exhibit 
a minimum. This is attributed to the crossover from the 
quartic-coupling dominant to the cubic-coupling domi- 
nant regime. As we will discuss later, the cubic term 
considerably stabilizes the glassy state and enlarges the 
region of glassy state in the phase diagram. This results 
in the initial drop of xiV values at the transitions as / 
deviates from 0.5 . 

For very asymmetric copolymers, mean- field theory 
predicts a first-order transition into ordered spherical 
phases [face-centered cubic (fee) or bed at smaller 
than the mean-field spinodal {xN)s However, ex- 

periments show that between the featureless disordered 
phase and the ordered bcc phase, there exists an inter- 
vening disordered-micelle state jH El IH IH HI El IH 
E3,EME3|- In a self-consistent-field calculation, it was 
shown |3(]| that the micelles are formed via a thermally 
activated process, with a free-energy barrier vanishing 
at the mean-field spinodal. Therefore, if the system is 
quickly quenched to the vicinity of the spinodal, micelles 
will proliferate all over the sample; the jamming of these 
micelles causes their translational diffusion to be so slow 
that long-range order cannot be developed. 

The interplay between the glass transition and the 
ODT is complicated. We will present some tentative re- 
sults in the next subsection. But here we simply note 
that, in contrast to the symmetric case where glass tran- 
sitions occur at > {x^)s (or below the mean-field 
spinodal temperature), for asymmetric copolymer glass 
transitions can occur at x^ < {x^)s (or above the 
mean- field spinodal temperature). We attribute this to 
the fact that different arrangements of micelles could gen- 
erate a large number of metastable states, which signifi- 
cantly stabilize the glassy state. 

As highlighted in Eq. the chain length N controls 
the magnitude of nonlinear fluctuations and, hence, the 
deviation from mean-field behavior that is recovered in 
the limit N ^ oo. Figure |21 shows the chain- length de- 
pendence of the glass transition temperatures [measured 
by {xN)a and {xN)k, respectively] relative to the mean- 
field spinodal for symmetric and asymmetric (/ = 0.3) 
copolymer melts. It is clear that in both cases the glass 
transitions (both Ta and Tx) approach the mean-field 
spinodal as N goes to infinity (though from different di- 
rections in symmetric and asymmetric cases), implying 
that in this limit the mean-field spinodal is true stability 
limit of the disordered phase (with respect to either or- 
dered or randomly phase-separated structures) . In other 
words, the mean-field spinodal is ultimately responsi- 
ble for the appearance of random structures, thus is the 
mean-field signature for the glass transition. This general 
conclusion is likely to be universal to the class of models 
with continuous degeneracy in the ground states, such as 
the Brazovskii model (see Ref. 66J for other models of 
the same class). We note that the connection between 
the spinodal and the glass transition was also implied 
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FIG. 2: Chain-length dependence of glass transitions. Plot of 
A (xN) = xN — {x^)s against iV. Dashed lines are for the 
Kauzmann temperature, and dash-dotted lines are for the on- 
set of glassiness. Upper dashed and dash-dotted lines are for 
the symmetric copolymer, and lower ones for the asymmetric 
copolymer with / = 0.3. 

in an earlier study by Bagchi and co-workers (b^ I of the 
Lennard-Jones liquid, where the authors conjectured that 
the liquid spinodal corresponds to the state of random 
close packing in an equivalent hard-sphere system. 

We now discuss the chain-length dependence of the 
gap between glass transitions and the spinodal, {x^)a — 
{xN)s and {xN)k — (x-^)s- By a simple scaling analysis 
given in Appendix lUl we find that both should scale as 
^-0.3 symmetric co poly mer. This is consistent with 
the analysis of Wu et al. if we substitute the N depen- 
dence of the parameters into their scaling relation. Our 
more accurate numerical calculations confirm this result, 
as shown in Fig. |3| where the first-order transition into 
lam phase is also included for comparison. 

For asymmetric copolymers, the results are more com- 
plicated: for short chains, both {x^)a and {xN)k are 
larger than the spinodal value {x^)s\ as N increases, the 
transition lines first shift downward below the spinodal 
line, which indicates a possible crossover; then, for even 
larger iV, the transitions gradually approach the spinodal 
and eventually collapse to the spinodal as iV — > cxd. In 
this latter limit, we find, using the scaling analysis out- 
lined in Appendix O that {xN)s - {xN)a,k '-^J^'^/^. 
We attribute the nonmonotonic dependence on N to the 
crossover from the quartic-coupling dominant to cubic- 
coupling dominant regime as the chain length increases. 
Generally for asymmetric copolymers, quartic coupling 
dominates for short chains and the glass-transition lines 
are located above the spinodal (or below the spinodal 
temperature); for long chains, the opposite holds. 

We close this subsection with a brief discussion of the 
dynamics of the system. For the Coulomb-frustrated- 
magnet model, by invoking the ent ropy-dr oplet picture, 
Wolynes and co-workers predicted [31, [T^ that the sys- 
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tern should exhibit relaxations similar to fragile liquids 
fflf , characterized by the Vogel-Fulcher behavior, with a 
relaxation time t cx exp[^/(r— Tr: )] between Ta and Tk, 
and diverging at the Kauzmann temperature Tr-. This 
prediction was disputed by Geissler and Reichman |65| . 
who performed dynamic Monte Carlo simulations of the 
Brazovskii model without the cubic interaction. They 
found that as the system approaches the glass-transition 
temperature predicted by the mode-coupling theory, the 
relaxation time indeed increases dramatically, but does 
not show characteristics of fragile liquids. Schmalian, 
Wu, and Wolynes subsequently argued [t^ that the fail- 
ure to find the expected dynamic behavior could be a 
result of the mode-coupling approximation which over- 
estimates the transition temperature. Here we note that 
the simulations by Geissler et al. were performed at tem- 
peratures above the ODT temperature, but our calcula- 
tions show that in the absence of the cubic interaction, 
the onset of glassiness always occurs below the ODT tem- 
perature, at least in the weak-coupling regime. Therefore 
simulations at lower temperatures (below the ODT tem- 
perature) are necessary in order to elucidate the dynamic 
behavior of this model. 

For block-copolymer melts, the situation is even more 
complicated. The Hamiltonian given by Eq. is a 
coarse-grained description of the system that focuses on 
the physics at length scales comparable to or larger than 
the size of the polymer. Therefore we expect the valid- 
ity of our analysis to be limited to this range of length 
scales. The configurational entropy Sc defined above only 
measures the number of configurations of chain aggre- 
grates in the locally phase-separated structures, but does 
not account for different chain conformations within each 
aggregrate. Indeed, above the glass-transition tempera- 
ture of the monomer, polymer chains remain liquidlike 
even though the system acquires solidlike behavior at 
the microstructural scale (at high frequencies). Chain 
diffusion also provides an additional mechanism for re- 
laxations. Therefore, to accurately describe the dynamic 
relaxations in block-copolymer melts, one has to consider 
relaxations both at the microstructural scale and of in- 
dividual chains. 



B. Glass-transition vs order-disorder transition 

Our analysis in the previous subsection shows that the 
glass transition is possible in diblock-copolymer melts 
and is related to the underlying mean-field spinodal of 
the disordered phase, which is responsible for the pro- 
liferation of inhomogeneous metastable states. However, 
the ODT also occurs in the neighborhood of the spinodal. 
Thus a full understanding of the glass transition in this 
system must address the relationship between these two 
transitions. 

In molecular liquids the glass transition always takes 
place in the supercooled state-below the melting (freez- 
ing) temperature. However, in diblock-copolymer melts. 
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FIG. 3: Glass transition vs ODT in symmetric copolymer 
melts. Dashed and dash-dotted lines have the same meanings 
as in Fig. |5| the solid line represents the ODT into the lam 
phase. 



the structural entities forming the random structures 
are themselves molecular aggregates formed through self- 
assembly, the number and size of which depend on the 
temperature of the system. Therefore the relationship 
between the glass transition and the ODT is not obvi- 
ous. 

Microphase transitions in block-copolymer systems 
have been extensively studied, both theoretically and ex- 
perimentally 0, 13 . In previous theories only the Hartree 
term arising from the quartic interaction was retained; 
fluctuations due to the cubic interaction were ignored. 
However, our analysis in the previous subsection shows 
that fluctuations due to the cubic interaction play an es- 
sential role in the glass transitions in asymmetric diblock 
copolymers, at least for long chains |73- Moreover, as 
discussed in Sec. Ill Bl the leading cubic diagram is of the 
same order in N as the Hartree term and their numer- 
ical magnitudes are comparable. Therefore we need to 
include fluctuations due to the cubic term in our studies 
to have a consistent comparison between the ODT and 
the glass transition. 

In this subsection, we compare the transitions into the 
ordered phase and into the glassy state. We have chosen 
to study symmetric (/ = 0.5) and highly asymmetric 
/ ~ 0.1 copolymers, as our perturbative methods are 
better controlled in these two limits (dominated by the 
quartic and cubic nonlinear interactions, respectively). 

Figure 13 shows the chain length dependence of the 
transitions for symmetric copolymers; it can be consid- 
ered as a generalized phase diagram. The solid line delin- 
eates the equilibrium phase boundary between the disor- 
dered phase and the ordered lam phase. For a given N, 
as temperature decreases [A{xN) increases], the equilib- 
rium state of the system will change from the disordered 
phase through a weakly first-order transition to the lam 
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phase. However, since the nucleation kinetics is gener- 
ally slow and complicated 111] , if the system is super- 
cooled to avoid the nucleation of lam phase, the system 
will remain in a metastable disordered state below the 
ODT temperature. Upon further cooling to the tem- 
perature Ta shown as the dash-dotted line, the system 
enters the glassy regime. The region bounded by this line 
and the Kauzmann line (the dashed line) defines the dy- 
namic range within which glass transition can take place 
[31I l3^ . Although the lam phase has the lowest free en- 
ergy at low temperatures, once a system is supercooled 
below Tfc, it becomes essentially frozen and incapable of 
reaching the more stable lam state. The narrow gap be- 
tween the onset of glassiness and the Kauzmann temper- 
ature implies that the glass transition in block copolymer 
melts will be fairly sharp. 

In symmetric copolymer melts we observe the scaling of 
{xN)oDT-{xN)s ~ as predicted in Ref. [l|. For 

the onset of glassiness, {xN)a — {x^)s scales as 7V~°-^, 
which agrees well with our approximate scaling analysis 
given in Appendix El Our results show that for sym- 
metric copolymers, the ODT always occurs before the 
glass transitions (i.e., at temperatures above the glass 
transitions). While one might argue that this conclusion 
could be due to the particular choice of diagrams in our 
perturbative calculation, we find that this scaling with 
TV remains unchanged when a different approximation 
scheme, the self-consistent-screening approximation, is 
used yu\ . In addition, our results are also consistent with 
the local-field calculations by Wu et al. ^Qj , as will be dis- 
cussed later in this subsection. Since {x^)odt ~ {x^)s 
decays more rapidly with N than both {xN)a — {xN)s 
and {xN)k — {x^)s, for sufficiently large N, we always 
have {xN)oDT < {x^)a,k- Therefore, at least in the 
long-chain limit, our conclusion that the glass transition 
occurs below the ODT temperature should be valid, re- 
gardless of the approximations in the calculation. 

Figure^shows various transitions for highly asymmet- 
ric copolymers around / = 0.1. Here again Ta is not 
shown as it is very close to Tk on the scale of the fig- 
ures. In the case ofN = 10^, the glass-transition lines are 
located below the ODT, i.e., the glass transition temper- 
atures are above the ODT temperature. In other words, 
glass transitions can precede the ordering transition into 
the bcc phase. This unusual behavior is quite different 
from what happens in molecular fluids, where the glass 
transition always occurs below the freezing (ordering) 
temperature. In the case of longer chains with N = 10*, 
the ODT occurs before the glass transitions; this is the 
expected behavior in the asymptotic limit N —^ 00, since 
in this limit the glass-transition lines approach the mean- 
field spinodal whereas the ODT into the bcc phase takes 
place at a finite distance below the spinodal |l3| . 

The chain-length dependence of the glass transitions 
relative to the ODT for asymmetric diblock copolymers 
is qualitatively similar to the critical micelle tempera- 
ture in the same system. It is shown [3(| that disordered 
micelles can appear in large numbers before the order- 



ing transition only for not-too-long chains; for very long 
chains, the ODT will set in before the disordered micelles 
reach a considerable concentration, essentially precluding 
the disordered micelles from being a distinct intervening 
phase between the featureless disordered state and the or- 
dered (fee or bcc) phases. Since micelles are likely to be 
the structural entities in the glassy asymmetric copoly- 
mer melts, the connection between the micelle formation 
and the glass transition is worth further investigation. 
As discussed in Sec. IIII Al the cubic interaction sta- 
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FIG. 4: Glass transition vs. ODT in very asymmetric copoly- 
mer, (a) iV — 10^, from above: mean-field spinodal (dotted 
line), ODT (solid line), and Kauzmann temperature (dashed 
line), (b) TV = 10*, from above: mean-field spinodal (dotted 
line), Kauzmann temperature (dashed line), and ODT (solid 
line) . 
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bilizes the glassy state. We attribute this stabihzing ef- 
fect to the additional complexities in the configurational 
space caused by the cubic term. This effect is closely re- 
lated to the effect of the cubic interaction on the ODT. 
Theoretical analysis shows that the presence of the cubic 
term can considerably reduce the free energy of ordered 
microstructures with three fold symmetries. This is con- 
sistent with the fact that there are more stable ordered 
phases in asymmetric diblock copolymers. If we visualize 
the random structures as polycrystals with local but no 
long-range order, then the increased variety of mesophase 
structures will increase the complexity in the configura- 
tion space [7^ , which can explain the stabilization of the 
glassy state in asymmetric copolymers. 

As a final technical point, we compare our treatment of 
the cubic term with that in Ref. (36|. There the authors 
used a local-field approximation, in which a momentum- 
independent self-energy is solved variationally by map- 
ping the Brazovskii Hamiltonian [as given by Eq. |(2J)] 
to a reference nonlinear but local Hamiltonian. Within 
this approximation, it was found that the cubic interac- 
tion considerably stabilizes the glassy state and the glass 
transition can occur at temperatures above the mean- 
field spinodal temperature; these results coincide with 
ours. However, for certain choices of parameters in the 
weak-coupling regime, this local-field treatment could re- 
sult in a nonmonotonic relation between the correlation 
length and the temperature. We believe this unphysical 
behavior is probably due to overestimating the fiuctua- 
tion effects due to the cubic term in their treatment. 



IV. CONCLUSIONS 

To conclude, using the thermodynamic replica formal- 
ism we have shown that at low temperatures, diblock- 
copolymer melts can exist as randomly microphase- 
separated structures, in addition to the thermodynam- 
ically stable periodic structures. This transition is essen- 
tially a glass transition in which the supercooled liquid 
gradually gets vitrified. We have identified the temper- 
ature range over which this glass transition can occur, 
which is bordered by the onset of glassiness (or the dy- 
namic glass transition) temperature from above and the 
Kauzmann (thermodynamic glass transition) tempera- 
ture from below. For symmetric diblock copolymers, the 
glass transition takes place below the temperature of the 
ODT into the lam phase. However, for asymmetric di- 
block copolymers, the glass transition can precede the or- 
dering transition, which is an unusual feature that prob- 
ably refiects the self-assembly nature of the system. This 
study leads us to naturally identify the quenched samples 
of block copolymers in some previous experimental works 



as the glassy state of the system. Given the slow phase 
transition kinetics in copolymer systems, we expect such 
glassy structures to be quite common in these systems 
without externally imposed aligning fields. 

As in any theories on polymer mixtures [l^ Isij , the 
scaled degree of polymerization, N, serves as a Ginzburg 
parameter which allows us to systematically examine the 
approach to mean-field behavior as N —^ 00. An im- 
portant conclusion is that in the limit of infinitely long 
chains, the glass transitions collapse to the mean-field 
spinodal, suggesting that the mean-field spinodal is ulti- 
mately responsible for the proliferation of inhomogeneous 
free-energy minima and can be used as the mean-field sig- 
nature for the glass transition. 

That a glass transition occurs at the mean-field spin- 
odal in the limit of ./V ^ 00 can also be understood using 
the following dynamical argument. Since the Hamilto- 
nian has an overall factor of 7V^/^, in the mean-field ap- 
proximation, we expect the free energy barriers between 
the metastable states to be proportional to this factor. 
For very long chains, these barriers can be very large. 
Since proliferation of the metastable minima appears at 
the spinodal [6^ . upon a quench below the spinodal, the 
system will first go to these metastable states with over- 
whelming probability because of their large number, and 
transitions from these metastable states should be very 
slow. Note that it is the barriers from these metastable 
states to the (more stable) ordered phases and between 
the metastable states themselves, rather than the nucle- 
ation barrier from the uniform disordered phase to the 
ordered phases, that are relevant to the glass transition. 
Hence, for example, in symmetric diblock copolymers, 
even though the transition from the disordered to lamel- 
lar phase approaches second order in the limit N ^ 00 
(where the nucleation barrier vanishes |46j|'). our theory 
predicts a glass transition that coincides with the ODT, 
which is the spinodal in this limit. 

Studying diblock-copolymer melt as a specific example 
of the Brazovskii model, we find that the cubic interac- 
tion significantly increases the stability of the glassy state 
as well as the bcc phase, and causes qualitative changes in 
the scaling relations with the chain length. We conjecture 
that this stabilizing effect is due to increased configura- 
tional complexity as a result of more free-energy minima 
due to the presence of the cubic term. 
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APPENDIX A: PERTURBATIVE EXPANSION OF THE EFFECTIVE POTENTIAL WITH BROKEN 

SYMMETRIES 

In this appendix we present the details of our perturbative calculation of the free energy defined in Eq. IjlTII . The 
general expansion of the effective potential for a system with broken symmetry was derived in Ref. "S^l . Here we omit 
the details of that derivation but give the result 

r[(^,G] =/[(p] + iTrlnG-i + iTr(D-iG) -r2[G;Ai/]. (Al) 

Here ip is the order parameter in the ordered phase, I[(p\ is the mean-field free energy (in our case the Leibler free 
energy), Ai/[^/;; if] is the shifted Hamiltonian [see Eq. 0], D is the shifted bare propagator defined as 

^ ^4^^, (A2) 

and G is the renormalized propagator. As noted before, we reserve boldface uppercase letters for matrices of correlation 
functions and use the corresponding plain ones when referring to the matrix element. The second term (Trln) in 
Eq. ljAl|l is the one-loop correction and the last term, r2 contains higher-order corrections, including all 2PI diagrams 
generated by the vertices in the shifted Hamiltonian Ai/ with the renormalized propagator G. The third term ensures 
the consistency of the expansion in terms of the renormalized propagator. 

It has been shown in Ref. [s^ that r[iy9, G] as defined in Eq. ljAl|l is stationary with respect to both ip and G. 
Therefore one can derive the self-energy equations through a variation of Eq. HA1|I . which gives 

In the field-theory description of diblock-copolymer melts [Eq. jSJ], iV~^/^ serves as a smallness parameter, which 
enables a straightforward loop expansion for r2 [G] . To the leading two-loop order (one-loop order in the self-energy) , 
one has three terms 

(A4a) 
(A4b) 
(A4c) 

corresponding to the diagrams shown in Figs. ISfa), 5(b), and 5(c) respectively. In the glassy state, translational 
symmetry breaking does not occur; therefore, = 0, Dabiq) — 9{<l)5ab^ and r2'^^ vanishes. Note that Fj^^ is the 

(2) 

Hartree term which only generates a momentum-independent self-energy in the diagonal part of G; F2 generates an 
off-diagonal self-energy which enables a nontrivial solution with broken replica symmetry. For symmetric copolymer, 
the cubic coupling is zero; therefore, to find possible solutions with broken replica symmetry we need to include the 
off-diagonal term of second order, 

^2^^ " 48^3/2 E / ^ '^^'>{qi)Gabiq2)Gab{q3)Gab{-qi - 92 - gs), (A4d) 

a.b 



r(i) _ 
2 ~ 


8N 


f d^qid^q2^ / / n 

1 j-2^-j6 '-'a.a[qi)l^aa[q2), 




p(2) 
2 — 




(27r)6 ^ab{qi)Gab{q2)Gabi 91 


92), 


p(3) 
^2 ^ 




r d^qi(Pq2<Pq-i , , , 

(27r)9 V'a(gi)Gafc( q2)Gab[ 


'q3)Gab{qi + 92 + q3)'fib{-qi), 



corresponding to the three-loop diagram as shown in Fig. ISJ^d). To study the crossover from very asymmetric to 
mmetric copolymer, we keep r2^-' in the off-diagonal 
From Eqs. HA3|I and (|A4|I we obtain the self-energy 



symmetric copolymer, we keep r2^-' in the off-diagonal renormalization for asymmetric copolymer as well 



Sa6(fc) = j -^Gaa{q)Sab " J -^^G ab{q)G ab{k ~ I) - ^ j ^^^Gab{q)Gab{p)Gab{k ~ P ~ q) ■ 

(A5) 

The three terms on the right-hand side corresponding to Fig.|SJe), 5(f), and 5(g), respectively. 
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FIG. 5: Feynman diagrams: (a)-(d) Loop diagrams in (e)-(g) Self-energy diagrams. In the diagrams thick solid lines 
represent the renormalized propagator G and wiggly lines represent the external leg of the order parameter (fi{q). We use a 
slightly different perturbative expansion for the diagonal renormalization, as explained in Appendix IHl 



Under the one-step-replica-symmetry-breaking (1-RSB) ansatz, Gab{q) ~ ^ ^{q}] Sab + ^{q), and r2 has three 
terms from Eqs. IIA4all . (|A4b|) . and l)A4d|l 



.(1) 

2 

(2) 
2 

(4) 



77 



mX 
2 



UN 
A2 



m 



i2n) 



6 G{qi)G{q2)G{~qi-q2) + m{m-l) I ^ ^^'^^'^ T(qi)T{q2)T{-qi - 92) 



48iV3/2 



(27r) 
Sqid^q2(Pq3, 
(27r)9 



(27r)e 



Q{qi)Q{q2)Q{q3.)Q{-qi - 92 - q-i) 



m{m-l) I ^ ^^(^2tt)^ •^(gi)-^(g2)-?^(g3)-^(-gi - g2 - gs) 



(A6a) 
, (A6b) 

(A6c) 



Using the polarization functions Iiab{k) defined as 



Tiabik) 

ng(fc) 
n^(fc) 



d3 



(27r)3 
(27r)3 



Gab{q)Gba{k + (?) 

- n^(fc)] Kb + 

g{q)Q{k + q), 



(2^)3 

[ng(fc)-n^(fc)]^,, + .F(fc), 



T{q)T{k + q). 



(A7a) 
(A7b) 



we can rewrite the self-energy functions (after taking m = 1) as 



Y.g{k) = 



{T,g — Sjr) Sab + Sjr, 

27V1/2 y (27r)3^^'^^ 



77^ /■ d3„ 



(A8a) 
(A8b) 
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The configurational entropy is obtained from Eqs. (|32() and |IA6|I to be 

^ 0(0) , c(l) 

= - ^ 



27V1/2 7 (27r)3 



7?^ 



A2 



d3q 



n^('z)n^(-'z)- 



UN J (27r)3 -"'^^ ' 48iV3/2 7 (27r)3 

We have set /csT = 1 in above derivations; the configurational entropy is given in unit of per unit volume. 



(A9a) 
(A9b) 



APPENDIX B: ORDER-DISORDER TRANSITION 



In this appendix, we present our calculation of the ODT in diblock-copolymer melts. Our approach is different 
from the Brazovskii approximation [isL l4ll| . 

Following the derivation in Appendix ^ we expand the effective potential to two-loop order and keep only the 
diagonal terms in Eqs. (jA4|) : 



1 



2iVi/2 



(27r)3 



Q{q) 



UN J (27r)^ 



g{q)Ug{-q)^ 



UN 



(Pqid?qi 
(27r)6 



V(9i)^(-92)ng(gi + q2)'p{-qi) 



(Bl) 



where FL[ip) is the Leibler free energy for the ordered phase, D{q) is the shifted bare propagator as given in Eq. ||UJ). 
The Hartree approximation (similar to the Brazovskii approximation) amounts to keeping only the first four terms 
of Eq. HB1|I . which can be justified by a renormalization-group argument 73]. The central idea is the following: 
since near the critical temperature the dominant fluctuations are those with wave numbers close to qm at which the 
propagator is maximized, one can decompose the spherical shell into small "patches" and rewrite the order parameter 
into n components, each corresponding to one patch. In this way one can rewrite the original Hamiltonian [Eq. |5J] as 
an n-vector model. At the critical point, n goes to infinity and the Hartree approximation becomes exact. Therefore 
we may replace Q by the Hartree approximation Qh as defined in Eq. (|12|l 

GHikr-Dik)-^^^,j^^g,iq). 

This gives the first three terms in Eq. H1U|) . 

However, here we want to study the correction due to the cubic coupling; thus, we want to include the leading-order 
diagram from the cubic interaction in the effective potential, as shown in Fig. lUb). It can be shown that in the 
corresponding n-vector model as mentioned above, the Hartree term [Fig. El^a)] is of order 0{n) and this correction 
term [Fig. El^b)] is of order O (n^/'^) . In our numerical calculations we find these two terms to be comparable for 
the temperature range we are interested in. By a similar argument the last term in Eq. IjBip is of order 0(1) and 
ignored in our calculation (the numerical value is indeed small compared with the other one-loop diagrams because of 
the weak first-order nature of the transition). To summarize, the free energy is given by Eq. I)B1|I with the last term 
dropped, as is Eq. fTT|l . 

To find the ODT temperature, the free energy is minimized numerically with respect to the magnitude of density 
wave a as given in Eq. Q and the ODT occurs when the free energy of the ordered phase equals the free energy of 
the disordered phase. 

The physical correlation function is given by Eq. (|13|l and the corresponding self-energy is 

^a(^) - ^ / ^.^Hiq) - ^ / ^,QHis)QH[q + k). (B2) 

This renormalization scheme includes two parts, the first corresponding to a simple Hartree approximation and the 
second incorporating fiuctuations from the cubic interaction using the Hartree-renormalized propagator, which is 
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schematically shown in the following diagrammatic equations: 




(B3) 

4=0- 



where thin lines represent the bare propagator g{q), double lines represent the Hartree-renormalized propagator 
Gh{q)^ and thick lines represent the physical propagator Q{q). Equation (|B2p modifies the self-energy equation (jASap 
we derived using a straighforward loop expansion in Appendix ^ One can verify that this self-energy equation does 
not have the unphysical nonmonotonic relation between temperature (manifested through in '''o) and correlation 
length [manifested in r in Q{q)\ which occurs in a naive loop expansion, and this renormalization scheme indeed 
gives consistent result in the known limits; e.g. when {xN)s — 3> 1 it reduces to the loop expansion and when 
(X-^)s ~ ~ it gives the leading-order terms in the 1/n expansion. 

APPENDIX C: APPROXIMATE SOLUTION OF THE GLASS TRANSITION 

In this last appendix we provide an approximate solution of the self-consistent equations obtained in Appendix fXl 
The diagonal and off-diagonal self-energy equations are shown in the following diagrammatic equations: 




where dashed lines represent the renormalized off- 
diagonal propagator JF; thick lines and double lines rep- 
resent the renormalized diagonal propagator Q and the 
Hartree-renormalized propagator, respectively, the same 
as before. 

From Eqs. ^ and lEH the renormalized propagators 
Q and J- are given by 



Gil) 



4g„ 



(28') 



(31') 



When r,s are small, the polarization functions can be 



approximated as 

ng(fc) 



1 

-(- 



(C3a) 
(C3b) 



for < |fc| < 2 and zero elsewhere. These are verified 
numerically and work well for r, s not too large (^ 0.1). 
The diagrammatic terms in our calculations are found to 
be 



_fq 



■ 5{q) 



27r^ 



{2n) 



-g{k-q)ngiq) 



k—Qr, 



Svrr-i 



(C4a) 



(C4b) 
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{2^ 



k—q-n 



-(--- 



(C4c) 



(C4d) 



— r3n^(-q)n^(q) ^ — 



(27r)3 



167r2 



. (C4e) 



From Eqs. (B3) and (28') we have the fohowmg equa- 
tions for r 



A 



To 



(C5) 



And from Eqs. (C2) and (31') we have 



A2 



487r7Vg2 



1 



Let us look at Eq. ljC6p first. Defining t 
Eq. I)C6|I becomes 



A2(l-<)3 



7?2(l-i)2 



487r7Vg2^r5/2 BA^ i/2g3^r2 ' 



1 

(C6) 

AT*: 



(C7) 



This equation always has a replica-symmetric solution 
t — 1 {T = 0). Here we are seeking a replica-symmetry- 
broken solution with t < 1. Defining the dimensionless 
parameters 



A 
B 



A2 



487r Ag^r5/2 ' 

r,2 



77' 



87Vi/2g3^r-2 ' 
Eq. HC7|I becomes 

A{l-t)H^ B{l-t)e 



l + t 



l + t 



1. 



(C7') 



Numerical calculations show that when both A and B 
are non-negative and either A > 23.66 or B > 11.09, 
there is always a solution < t* < 1. For symmetric 
copolymer and very asymmetric copolymer, respectively, 
these inequalities result in the criteria 



< 



< 



A2 



2/5 



7V-2/5, 



1/2 



-1/4 



(C8) 
(C9) 



And the resulted scaling relations for tq [oc (x-^)s — X^] 
for symmetric and asymmetric copolymers are, respec- 
tively, 

-N-°-^, (CIO) 
Af-i/4. (CU) 



To 
To 



These have been verified by our numerical calculations. 

Finally we look at the configurational entropy and the 
Kauzmann temperature. The configurational entropy is 
given in Eqs. HA9|I and found to be 



r(l-0= 



(C12a) 



967rA^r3/2 ^ ' YQg^2]y3/2^2 

(C12b) 

Thus for symmetric copolymer (77 = 0), the Kauzmann 
transition is located at 

rr^N-^/^, (C13) 
To - (CM) 

And for very asymmetric copolymer {rj/qm^ 3> X/qm), 
r^N-^/^, (C15) 
To ~ (C16) 



APPENDIX D: RELATIONSHIP BETWEEN THE 
PINNED FREE ENERGY F[(] AND THE 
FREE-ENERGY LANDSCAPE OF THE 
ORIGINAL HAMILTONIAN H[(l>] 

Here we explicitly show that the free energy F[C,] de- 
fined in Eq. H14|l captures the metastable free-energy min- 
ima of the Hamiltonian H[(j)\ as defined in Eq. (|2Jl. First 
we rewrite Eq. (fTH) as 



where * is a shorthand notation for integration and 

H'[^]=H[^] + ^<j>*<j>. (Dl) 

We then define the generating functional of the perturbed 
Hamiltonian H'[(j)], 

W[J = aC] = - In J V(f>exp (-7^'[0] + J * 0) . (D2) 

The effective potential T'[(p] of the Hamiltonian H'[(j)] is 
obtained as the Legendre transform of VF[J]. Thus we 
have 



SW[J] 



SJ 



J=aC 

r'M 



5T'[ 



6(p 



= -V. (D3) 
= VF[J] + J*V5, (D4) 
= J. (D5) 
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Now W[J] is related to F[C] by 

W[J = aC] = F[C] - f C * C, 
so that for any (.* that minimizes F[Q, we have 



that is, 



SW[J] 



SJ 



+ aC = 0, 



J=qC* 



Equation ljD8|) holds for any positive a, including in par- 
ticular the limit a ^ 0+ . In the limit of a ^ , H' [(f)] 
(D6) approaches H[(j)] and r'[ip] approaches F [(/?], the effec- 
tive potential of the original Hamiltonian -ff [</)]. Also 
J = Qf^ ^ 0, from Eq. (|D5|I . (p* becomes a minimum 
of T[ip]. This, together with Eq. (|D8|I . shows that the 
minima of coincide with the minima of the effective 
potential T[ip] of the orginal Hamiltonian in the limit 
a 0+. 



6W[J] 
SJ 



if . 



(D8) 
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